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We study the Casimir force in a film with free boundary conditions described by the n -> oo limit 
of the 0{n) models. Via exact calculations we reduce the model to an one dimensional model the 
properties of which we determine numerically, near the critical region, and in an exact analytical 
manner in the low-temperature regime. The results resemble the experimental curve behavior for 
the force obtained for ''He films. 
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I. INTRODUCTION 

A recent paper by Diehl et. al. [T] reports on a study 
of the Casimir force in thin films (i.e. dimensions co^ x L) 
of a Ginzburg-Landau- Wilson (GLW) version of the 0{n) 
model in the limit n — >■ cxd, the system being subject to 
free boundary conditions in the finite direction. Because 
of the continuous symmetry of the model, which is bro- 
ken at low temperature when L — > oo, as well as the fact 
that the boundary conditions correspond to those that 
are appropriate in the case of ''He films the results of 
the calculations in [1, are relevant to the Casimir force 
measurements on such films described in [21 [3] . The su- 
perfluid transition in *He is, of course, correctly modeled 
in terms of the XY , or 0(2), model, and the results in 
refs. [21 [3] have been quite successfully reproduced by 
Monte Carlo simulations of this model in [4j and |5j . 

Nevertheless, the 0(n — oo) model merits consider- 
ation as a depiction of systems with broken continuous 
symmetry in the bulk insofar as this model is susceptible 
to a combination of analytical and numerical approaches, 
yielding both quantitative and qualitative insights into 
the behavior of those systems. In the current article we 
study the n — > oo limit of the 0{n) models without us- 
ing the mapping on the GLW model, but relying on the 
definition of the corresponding spherical model [6j; we 
recall that the infinite translational invariant standard 
spherical model is equivalent to the n — >■ oo limit of the 
corresponding system of n-component vectors [7l [8] , but 
for the spherical model with surfaces or, more generally, 
without translation-invariant symmetry, this equivalence 
is preserved only if one imposes spherical constraints in 
a way which ensures that the mean square value of each 
spin of the system is the same [S] — that is, one averages 
thermally, but not spatially. 

Generally such a model is considered analytically in- 
tractable. However, as we demonstrate here, this model 
can be analytically reduced via exact calculations to an 
one dimensional model the properties of which can be 
then either studied numerically, near the critical region, 
or in an exact analytical manner in the low-temperature 
regime. The differences of our approach from the one 
presented in Ref. [1] are i) the model is self-contained 
and does not require the determination of any parame- 



ters, e.g. the strength of the fourth order coupling in the 
GLW model, and ii) there is no need to use an effective 
width of the system, keeping in mind that there is always 
some small (microscopic) ambiguity in the definition of 
the thickness L; see, e.g., Ref. [TU] and the discussion 
therein. Finally, the present formulation of the model is, 
unlike the GLW approach, suitable for investigation of 
the properties of the system at all temperatures, particu- 
larly at temperatures considerably below that of the bulk 
transition and thus well outside of the critical regime. 

It is worth reiterating that the Casimir effect remains 
the object of intense studies, in its original formulation 
due to Casimir [11] (see the reviews [HI [13]), and in its 
thermodynamic manifestation [14j — see, e.g, the general 
reviews [151 |TB] and the reviews devoted to some specific 
aspects of the critical Casimir force [TT] [TH] . 



II. DEFINITION OF THE MODEL 

For an 0(n),n > 1 model of a d-dimensional system 
at a temperature T and geometry oo''"^ x L the thermo- 
dynamic Casimir force is defined by [T^], [TS] 



^Casimir t^) ^ qj^ ' 



(2.1) 



where fil\T,L) is the excess free energy 



/W(T,L) = /(^)(r,L) - L/fc(r), (2.2) 

and the superscript r denotes the boundary conditions. 
Here f^'^\T,L) is the full free energy per unit area of 
such a system subjected to the boundary conditions t 
and is the bulk free energy density. 

Consider a c?-dimensional square lattice, each lattice 
site occupied by an n-component classical vector spin 
having ferromagnetic interactions with its nearest neigh- 
bors. We single out one dimension, z, to be L lattice 
spacings long. At each of the L sites along the finite di- 
mension, there is a (d — l)-dimensional transverse layer 
containing a total of A spins, where A is large and will 
later be taken to infinity. Free boundary conditions are 
imposed in the z direction by putting a layer of zero 
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length spins on the top and the bottom of the film (i.e. 
at z = and z = L + 1). 

The model, as described, is not especially amenable 
to analysis. However, in the n — ^ cx) limit it is 
equivalent [3 [3] to a form of the spherical model, wherein 
the vector spins are replaced by real-valued scalar spins 
and each {d — l)-dimensional layer satisfies an individ- 
ual spherical constraint ^ = A. We simplify matters 
further by using the mean spherical model, which gives 
the same results as the spherical model in the thermody- 
namic limit, i.e. with A ^ oo. 

Our Hamiltonian is therefore 

H = -jJ2ss' + jJ2AAY.'h^A (2.3) 

where the first summation is taken over nearest neighbor 
spins s and s' , which may each lie in the same layer or 
in adjacent layers. J > is the ferromagnetic coupling 
and Ai is the spherical field for layer i which will be used 
to enforce the mean spherical constraints. The notation 
Si j refers to spin j in layer i, with i — 1,...,L and 
j = l,...,A 

Fourier transforming spins according to 



(2.4) 



where the sum extends over the first Brillouin zone of 
layer i, we find 

H = -JAJ2^^ + ^^s(q)t?^(q)s(q) (2.5) 



with 



and 



d-l 



cos Qk (2.6) 



k=l 



M,, 



2A„ i=j 

-1, N-Jl = l 
0, otherwise 



(2.7) 



After computing the partition function in the standard 
way, we find that the free energy per transverse unit area, 
in units of fc^T, to be 



~A 



+ ^5]ln[det(7^(q))]. (2., 



The spherical constraint is enforced in the mean via the 
Lagrange multipliers A^. In particular, we must have 







for each i — 1, . . . , L. 



III. RESULTS ON THE MODEL IN d = 3 

We specialize now to the case of three dimensions. We 
write using an outer-product expansion. Taking 

{a/,|V'^'^)} to be the eigenvalues and normalized eigen- 
vectors, respectively, of the matrix defined in Eqn. 
(2.7), we have 



H(q) 



E 

1=1 



a; — 2 cos Qj. — 2 cos qy 



(3.1) 



In the transverse thermodynamic limit, A — )• oo, the 
sum over Brillouin zone is reproduced by an in tegral , i.e 



I^q ^ (2¥F / '^^'i- ^^Sht of that, Eqns. (2.9 1 and (3.1) 



give three-dimensional spherical constraints 

L 



1=1 



ai 



ai 



(3.2) 



for each i = 1, . . . , L, where K{k) is the complete elliptic 
integral of first kind with modulus k. 

Furthermore, the free energy can be written in closed 
form for c? = 3: 



A 2 



(3J 
2n 



1 ^ 



1=1 L 



In a; - 2/3 JA, 



11 2 2 2-1^ 

I 



2' 2' 



at 



(3.3) 



with generalized hypergeometric function The bulk 
free energy per length, /;,, is also computed to be 



/3/b = -13 J A + I In 



2tt j 



1 

J dq^: dqy sinh"^ ■ 



(3.4) 



with u — i ■y4^^Tcos~^^7^^T<XDsq^ and A = 3 for /3J > 
PcJ, while A for /3J < jS^J is determined by the bulk 
constraint 



/3J = 



d^c 



A — cos q^ — COS qy — cos q^ 



(3.5) 



where bulk critical coupling /3c J is provided by the above 
equation via setting A = 3 in it. 

Therefore, given a temperature, the problem is to solve 
Eqs. (3.2) simultaneously, with i = 1,...,L, for the 



unknown spherical fields. A.;. Then, we construct the 
Casimir pressure 



A 



dL [ A 



PF{L + 1) /3J"(L - 1) 



A 



A 



(3.6) 
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Implementing this analysis numerically, we find the 
Casimir force curves shown in Fig. [l] for several sys- 
tem sizes L. The scaling variable \s x — {L/^^ jt, 
where t = {T — Tc)/Tc is the reduced temperature and 
^f,(i — > 0+) — £,Qt^'^ is the bulk correlation length. Here, 
I' = 1 is the corresponding critical exponent in the bulk 
spherical model in d = 3, and = [47r(/3j, the 
non-universal amplitude as determined from earlier re- 
sults [inj. As L — 7^ cx), the scaling corrections van- 
ish and we are left with the Casimir scaling function 
■d{x) = L^l3Fca.s{x)/A for this system under free bound- 
ary conditions. In fact, the curves for L = 200 and 
L = 500 are substantially the same, so the solid black 
L — 500 curve is already an excellent approximation to 

We note here that the full range of negative x appro- 
priate to each curve depends on L; the T = limit is 
achieved at x = —L/^^. 

From our numerical results, we find a Casimir ampli- 
tude of 



A 



-0.0107737 



and an extremum of the scaling function 



-0.1269001. 



(3.7) 



(3.8) 



These quantities are expected to be universal and, in- 
deed, they agree to great precision with the values given 
by Diehl, et. al. [T] for their closely related model. 
The location of the Casimir force extremum is at a:: = 
X = —8.08 for L = 500. Agreement with the measure- 
ments in [21 13] is less satisfactory, as it should be given 




FIG. 1. (Color online) The scaled Casimir force L^^FcasM 
as a function of scaling variable x = {L/^Q)t, with reduced 
temperature t = {T — Tc)/Tc and bulk correlation length am- 
plitude ^q", for L = 10 (top, solid blue), L — SO (second 
from top, dashed purple), L = 200 (second from bottom, 
dot-dashed red), L — 500 (bottom, solid black). The zero- 
temperature value of —({3)/8n is indicated as a horizontal 
dashed line. The L = 200 and L = 500 curves lie on top 
of eachother in the critical region a:: ~ 0, and both are close 
approximations to the scaling function '&{x). 



the widely-acknowledged difference between the 0{2) and 
0(oo) models. 



IV. ASYMPTOTIC LOW TEMPERATURE 
RESULTS IN d = 3 

The system may be studied in closed form for T ^ Tc, 
i.e. (3J ^ f3cJ. In that regime, the left hand side of 
Eqn. (3.2) becomes large, forcing the lowest eigenvalue. 



fli, to approach 4 and dominate the right hand side. This 
effectively renders the constraint as 



[3 J 



1 



27r fli V fli 



2 4 



(4.1) 



(1) 

As this must hold for each i, we will have il^] ' — const. 
Numerically, we see that this eigenvector approaches a 
constant vector (as opposed to one which varies with 
respect to i). In order for the constant vector to be 
an eigenvector of Ai with eigenvalue 4, we must have 
A = A* = (5, 6, 6, . . . , 6, 6, 5)/2. We now expand about 
A = A* using perturbation theory. 

The eigenvalues and normalized eigenvectors of A4 (A* ) 
may be found in closed form as 



A; = 4 -I- 4 sin^ 



7^(^-1) 
2L 



(4.2) 



and 



a(0 



cos I -^^]^(2i - 1) ) /norm, (4.3) 



where I = 1, . . . , L and the norm is y/L for I = 1 and 
\J Ljl for all other I. Applying perturbation theory and 
using the constraint, Eqn. (3.2), we find 



ai « 4 32e 



-47rL(fl-fle) 



(4.4) 



while 



a; 



27rL(i? 



i^(4/A^+2-0 



0[R- 



(4.5) 



for I > 2, where we now abbreviate R = l3J and Rc 
the corresponding value at the bulk critical point. As 
before, K{k) refers to the elliptic integral with modulus 
fc, and A„ is given by Eqn. ( |4.2[ ) even when n is not 
in the range 1, . . . ,L. Note that there is one eigenvalue, 
fli, which exhibits markedly different behavior from the 
other eigenvalues. 

In order to compute the Casimir force, we must de- 
termine the effects of th e perturbed eigenvalues on the 
free energy, Eqns. (3.3 1 and (3.4). A careful analysis. 



whose details are deferred to a later publication, shows 
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FIG. 2. (Color online) The scaled Casimir force (symbols) as 
compared to the closed form asymptotic results (solid curves) 
presented in section |IV] for L = 50, 100, 200 and 500. The 
asymptotic results turn out to be accurate for moderately 
low absolute temperatures, corresponding to reduced temper- 
atures t < —0.2 in this plot. Note that the plot is versus un- 
sealed reduced temperature, t, which is equal to -1 at T = 0. 



that there is a leading order, temperature-mdependent, 
term in /3i^cas- This term is given by 



1 



where 



167r2 



sinh ^ 



d\ hi 



1 - £-4(^+1)1' 
1 - 6-4(^-1)-" 



-y/A — 2 COS Qx — 2 cos Qy 



(4.6) 



(4.7) 



Expanding in powers of 1/L, we find 



c as 



87rL3 



C(3) + 



2C(3) + iC(5) 



+ 0(L-) 



+ 0{T). (4., 



The leading temperature-dependent term is of the form 



4(i? - R,)L^ 



a + b\nL + 0(L- 



(4.9) 



with 



C'(-2) 



3K 



211n2-|-61n(27r) 
3C"(-2) 



(4.10) 



and 



b = 



3C'(-2) 



(4.11) 



where C is the Riemann function. Note that l/(i? — 
Rc) T for r « 0. The comparison between numeri- 
cal and asymptotic results is shown in Fig. [2j using in 
addition the next order results in powers of L which are 
omitted here for brevity. 

While the Casimir forces in Figs. [T] and j2| differ in de- 
tail from the data for ^He obtained in [21 13], the overall 
features — particularly the pronounced minimum in the 
Casimir force below the critical point and the approach 
to a non-zero limiting value at low temperatures — are 
strikingly similar. This encourages us to anticipate that 
the model will prove to be a very useful adjunct to gen- 
eral, and perhaps specific, studies of Casimir forces in 
systems with a continuous symmetry of the type that is 
broken in the superfluid transition. 
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